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A combined BATSE-GBM hard X-ray catalog is presented based on Earth Occultation Imaging applied to a 
reanalysis of BATSE data. An imaging approach has been developed for the reanalysis of Earth Occultation 
analysis of BATSE data. The standard occultation analysis depends on a predetermined catalog of potential 
sources, so that a real source not present in the catalog may induce systematic errors when source counts 
associated with an uncatalogued source are incorrectly attributed to catalog sources. The goal of the imaging 
analysis is to find a complete set of hard X-ray sources, including sources not in the original BATSE occultation 
catalog. Using the imaging technique, we have identified 15 known sources and 17 unidentified sources and 
added them to the BATSE occultation catalog. The resulting expanded BATSE catalog of sources observed 
during 1991-2000 is compared to the ongoing GBM survey. 



1. Introduction 

The past 20 years have seen unprecedented growth 
in the field of gamma-ray astronomy. Highly success- 
ful missions such as the Compton Gamma-Ray Ob- 
servatory (CGRO) detected gamma-ray sources and a 
more fundamental understanding of the basic physi- 
cal processes involved for those sources. From 1991 
to 1995, the Burst and Transient Source Experiment 
(BATSE) on CGRO, which was sensitive from 23 keV 
to 1.8 MeV, was the only instrument which was ca- 
pable of monitoring the low energy gamma-ray sky. 
In 1995, the Rossi X-Ray Timing Explorer (RXTE) 
satellite was launched, but it had a relatively narrow 
field of view and was only sensitive up to ~ 250 keV. 
Gamma-ray Burst Monitor (GBM) currently flying on 
the Fermi spacecraft since 2008 is the only instrument 
performing all-sky observations above ~ 200 keV and 
below 20 MeV at present. 

Therefore, the combination of BATSE and GBM 
all-sky analysis gives us a ~ 20 year history of the 
hard X-ray /soft gamma-ray sky. 



2. Earth Occultation Technique 

As demonstrated with BATSE [I] [2] and GBM [6], 
monitoring known hard x-ray /soft gamma-ray sources 
is possible with simple non-pixellated or low spatial 
resolution detectors by using the Earth as a temporal 
modulator of the gamma-ray flux. This is the Earth 
Occultation Technique (EOT). The Earth, seen from 
a BATSE detector, has a diameter of ~ 140°. When 
the source sets below the Earth's limb, a downward 
step appears in the detector count rate (and vice- versa 



for rising steps). The height above the background of 
the steps gives a direct measurement of the count rate 
of the occulted source. 

Occultation times can be predicted using the coor- 
dinates of the source and the spacecraft position. The 
occultation time is defined to be the time when the 
transmission through the atmosphere is 50% at 100 
keV. 

The Enhanced BATSE Occultation Package 
(EBOP) was developed by the Jet Propulsion Lab- 
oratory (JPL) group, and was applied to BATSE 
data to search and fit the occultation steps [2 . For 
known sources, a fitting procedure can be used on a 
240 second fitting window centered at the expected 
occultation time bin. For each detector and energy 
channel the model fit consists of a quadratic curve 
describing the background and a set of response 
vectors that model the effects of all sources contained 
within the fit window. However, the sensitivity of the 
step searching technique is limited by the difficulty 
in accurately determining the background over 
several orbits and by non-statistical noise introduced 
by bright pulsating sources with periods of a few 
minutes, such as Vela X-l and the unknown sources. 
Hence, unknown sources affect the EBOP results 
systematically. 



3. Earth Occultation Imaging Technique 
with Differential Filter 

Since EOT and EBOP require an input catalog 
with predetermined source position, another method 
is needed to search for unknown sources. An imag- 
ing method conceptually similar to x-ray computed 
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tomography used in medical imaging has been devel- 
oped 

When a source occults, the Earth's limb can be pro- 
jected onto the sky over which counts are integrated. 
During the orbital precession of the spacecraft, the an- 
gle between the source position and the plane of the 
satellite orbit changes. Thus over an orbital proces- 
sion period a range of projections are sampled. The 
projections from rising and setting occult ation steps 
can be combined to generate an image of the sky and 
localize a source. 

The angular resolution of this method is limited by 
the finite time required to go from full transmission 
to 50% attenuation for a source. For a typical or- 
bital speed of 8km/ s, the duration of an occultation 
step is ~ 8 seconds for normal incidence. As the 
angle between the orbital plane and the source in- 
creases, the duration of the occultation step increases 
by ~ 1/ cos(/3), where f3 is the elevation angle be- 
tween the source and the orbital plane of the satellite. 
The angular resolution is 360° x At/P, where At is 
the occultation duration and P is the satellite's or- 
bital period ( ~ 90 minutes). The angular resolution 
ranges from - 0.5° at f3 = 0° to - 1.25° at /3 = 66°. 
When /3 > 66°, occultation no longer occur. 

Instead of a catalog of known source positions, 
our imaging technique starts with a catalog of vir- 
tual sources covering the whole sky with spacing of 
~ 0.25°, resulting in a catalog of ~ 660,000 positions 
on the sky. With these virtual source positions and 
the spacecraft position history, occultation times can 
be predicted. To be able to combine rise and set steps 
together, occultation windows for rise steps are con- 
verted to set steps by rotating the data about the 
occultation time at the centre of the 4-min window. 
A weighted average of all the windows for a virtual 




Figure 1: Top: Example of an occultation window for the 
Crab. Bottom: The filtered data with the central portion 
fit to a polynomial (green), and the outer background 
portion fit to a spline function (red). 



source is calculated for an orbital procession period. 
Then, the averaged window is passed through a dif- 
ferential filter of the form [5], 

. = i^±fa ^ U=i-U \ (1) 

fb 

where Tj is the number of counts in bin j, f a is 
the inner boundary of the filter, and /& is the outer 
boundary (citation needed). The values of f a and 
used here are 3 bins and 8 bins, respectively. The dif- 
ferential filter is an approach to transform those step 
features into gaussian-like peaks, where the amplitude 
of the peak is related to the intensity of the source. To 
calculate the amplitude, the window is fit in two sec- 
tions. The central part, within ±(2/ a + fb) bins from 
the middle, is fit to quadratic, cubic, and quart ic func- 
tion, and the fit giving the lowest reduced x 2 will be 
kept. The outer portions of the window are fit with 
a spline function and joined by a straight line cross- 
ing the central region. The amplitude of the virtual 
source is determined by taking the difference between 
the two fits at the occultation time. 

Fig. [I] shows an example of an occultation window 
and the result of filtering for a virtual source very 
near the Crab. In fig. [I] the middle dotted vertical 
line corresponds to the expected occultation time of 
the virtual source. 



4. Data Selection and Source Confusion 

Prior to application of the Earth occultation 
method, proper data selection is required to remove 
large fluctuations that may affect the imaging. This 
selection process contains three steps. The first step 
is performed by using Quality Flags. Any high-rate 
events, which triggered the GRB acquisition mode, 
are flagged. The second step is due to South Atlantic 
Analomy (SAA) passage, where data bins within 20 
seconds before and after SAA passage are excluded. 
The last step consists of additional flagging of very 
short cosmic-ray events, which normally occur in only 
one detector at a time. 

During the Earth occultation imaging process, an- 
other systematic effect is source confusion caused 
by bright sources. The issue of source confusion 
arises from the location ambiguity of where along 
the Earth's limb the measured amplitude originates. 
Since BATSE detectors do not have imaging capa- 
bility, the flux for an occultation time could be al- 
located to anywhere along the Earth's limb. Hence, 
bright sources cause "X" patterns that can extend for 
tens of degrees. These "arms" make it difficult to 
search for fainter sources along the "arms". To mini- 
mize the effects of source confusion, an algorithm has 
been developed to ignore occultation windows for vir- 
tual sources when bright sources occult close in time 
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Figure 2: Example of Source Confusion "Arms". Source confusion algorithm has been applied to Crab and Cyg X-l, 
but not Sco X-l nor Vela X-l. Obvious pattern differences are shown in the map. 



(< 11 seconds). Fig. [2] illustrates the result of source 
confusion. The algorithm is applied here to the Crab 
and Cyg X-l. Hence, no bright arms originated from 
these two bright sources. But Vela X-l and Sco X-l 
are associated with the strong "X" patterns, since the 
algorithm was not applied here to them. 



5. Results 

The imaging algorithm has been applied to 64 pre- 
cession periods, 9 years of data, from 1991 to 2000, 
to generate all-sky images at 23-98 keV. Fig. [5] is 
a significance map with green contours at 3 a. In 
Fig. [5} the source names in black denote the origi- 
nal sources in the 2007 catalog; blue denotes the 32 
new sources; magenta marks 11 blank sources included 
to study systematic effects and sensitivity. These 
sources were identified through cross-correlating fea- 
tures in the maps with sources in the Swift/BAT, IN- 
TEGRAL/SPI, and Fermi/GBM catalogs. 

126 sources are included in the new BATSE/EBOP 
2012 catalog. Of those, 94 sources coincide with 
sources in the 2007 EBOP catalog, and 32 are new 
sources, with 15 of them identified with known Swift 
and INTEGRAL sources, and 17 unidentified. Ta- 
ble. [I] shows the name and position of the new identi- 



Table I 15 New known Sources found with Earth 
Occultation Imaging. 



Source Name 


RA 


Dec 


MCG+08-03-018 


20.67 


50.07 


RBS0315 


36.26 


18.81 


MRK595 


40.46 


7.19 


MRK3 


93.94 


71.03 


NGC3227 


155.87 


19.86 


NGC4593 


189.91 


-5.35 


NGC5506 


213.31 


-3.21 


PKS1510-08 


228.22 


-9.09 


IGRJ18485-0047 


282.17 


-0.79 


4U1850-086 


283.27 


-8.71 


MRK509 


311.03 


-10.73 


SWIFTJ2118.9+3336 


319.86 


33.56 


NGC7172 


330.51 


-31.87 


4U2206+54 


331.96 


54.51 


NGC7582 


349.60 


-42.37 



fied sources. 

Comparison with the current GBM catalog [6], it 
shows 69 of the BATSE sources coinciding with the 
GBM catalog. Of these, 3 are new sources found by 
the current imaging analysis, NGC5508, PKS 1510-08, 
and 4U 2206+54. 
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Table II 17 Unidentified BATSE Sources found with 
Earth Occultation Imaging. 



Unidentified Source 


RA 


Dec 


EBOP061-11 


61.19 


-11.70 


EBOP063+62 


63.37 


62.55 


EBOP100-9 


100.24 


-9.17 


EBOP168+18 


168.40 


18.40 


EBOP248-33 


248.39 


-33.92 


EBOP252-01 


252.02 


-2.19 


EBOP258-13 


258.30 


-13.26 


EBOP261-6 


261.63 


-6.93 


EBOP264-19 


264.47 


-19.53 


EBOP275+26 


275.37 


26.49 


EBOP281+55 


281.38 


55.51 


EBOP281-16 


281.41 


-16.05 


EBOP291-15 


290.82 


-15.01 


EBOP294+10 


294.25 


10.01 


EBOP305+6 


305.23 


6.00 


EBOP306+20 


306.03 


20.29 


EBOP346-8 


346.18 


-8.69 
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Figure 3: All-sky map summed over 64 precession periods (TJD 08393 - TJD 11690) over the energy range 23 - 98 
keV. 3 sigma (statistical) features are contoured in green, he source names in black denote the original sources in 
the 2007 catalog; blue denotes the 32 new sources; magenta marks 11 blank sources included to study systematic 
effects and sensitivity 
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